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Spherical symmetry for /(i?)-gravity is discussed by searching for Noether symmetries. The 
method consists in selecting conserved quantities in form of currents that reduce dynamics of f(R)- 
models compatible with symmetries. In this way we get a general method to obtain constants of 
motion without setting a priori the form of f(R). In this sense, the Noether symmetry results a 
' physical criterium. Relevant cases are discussed. 
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o ■ 

£h ■ I. INTRODUCTION 

The issue to renormalize General Relativity (GR) leads to effective theories of gravity including corrections in the 
curvature invariants to the Einstein-Hilbert action [lHj] ■ In addition, the discovery that the expansion of the universe 
is accelerated [5J has led to the introduction of the so-called dark energy to explain the accelerated expansion of the 
Hubble fluid which, up to now, has no agreed interpretation at fundamental level. 

In order to avoid this ad hoc ingredient, it has been proposed @ that, perhaps, one needs to modify GR at the 
O" 1 cosmic scales to address the observational evidences. 

5^ . In particular, /(i?)-gravity, where / is a function of the Ricci scalar R (in particular an analytic function), has been 
considered to explain the present cosmic acceleration without dark energy (see 

for reviews). Here we focus on 
the metric version of these theories, where the field equations are of fourth order for the metric. 
I ' The full investigation of such theory is the first step toward the so-called Extended Theories of Gravity where 
J> , effective actions involving generic curvature invariants are taken into account The goal is to address, by 

gravitational sector, problems like dark energy and dark matter at infrared scales and renormalization at ultraviolet 
| scales. 

In this view, exact solutions of the field equations are important in order to gain insight into the mathematical 
and physical content of a theory [l^|. They can be obtained only upon assuming special symmetries. In the case 
of metric /(i?)-gravity in vacuo, assuming spherical symmetry does not automatically lead to the Schwarzschild or 
Schwarzschild-de Sitter solutions [17j . The latter do solve the f(R) vacuum field equations, but the Birkhoff theorem 
of GR does not hold in metric f(R) gravity (see [Tel , [l9j for a detailed discussion) and, therefore, these are not all 
the solutions within the family of spherically symmetric solutions (20l - [23| . In fact, the class of spherically symmetric 
solutions of metric /(i?)-gravity is still unexplored, with only a handful of analytical solutions beyond Schwarzschild- 
de Sitter being known. When one considers axially symmetric solutions, the situation is even worse: apart from the 
Kerr metric [24j, only the solutions of Ref. |25fl are presently available. 

Because of the fact that these models have fourth order field equations, it is usually difficult to find exact solutions 
in full generality without imposing a-priori a particular form of the f(R) function. To this aim, it is possible to 
consider spherically s ym metric background and demonstrate that it is possible to find exact solutions via the Noether 
Symmetry Approach |l9l. [26rl28| . In fact, choosing an appropriate /(i?)-Lagrangian, it is possible to find out conserved 
Noether currents which will be useful to solve dynamics. This approach is very powerful due to the fact that it allows 
to find a closed system of equations, where we do not need to impose the particular form of the f(R) which, on the 
other hand, is selected by the Noether symmetry itself. In Section [IIJ we introduce the /(i?)-gravity using a suitable 
Lagrangian approach in the hypothesis of a spherically symmetric metric. Dynamical equations are solved in the case 
of constant curvature recovering the standard case of the Schwarzschild-de Sitter solution. In Section Hill we describe, 
in general, the Noether Symmetry Approach and in Section lIVl we apply it to spherically symmetric /(i?)-gravity. 
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In Section [V] we find solutions to the dynamics for both constant curvature and in the hypothesis of a static metric. 
Conclusions are given in Section [VTl 



II. SPHERICALLY SYMMETRY IN f(R) GRAVITY 
Let us start with a generalization of the Hilbert-Einstein action 



S 



d 4 x^f(R), 



(1) 



where f(R) is an analytic function of the Ricci scalar R and g is the determinant of the metric g^ v . We consider a 
time dependent spherically symmetric metric 



ds 2 = e^dt 2 - e x ^dr 2 - rW 



(2) 



where f(r,t) and A(r, t) are unknown functions of r and t and dil 2 — (dO 2 + sin 2 8dcf) 2 ), using unit of c = 1. In the 
following equations, we omit the t and r dependence, to avoid burdening the text. The Ricci scalar can be written as 



R 



1 2 

v rr H — v r (y r — A r ) H — {y r — X r ) 



-\ t {v t - A t ) - At 



(3) 



where the subscripts r and t indicate radial and time partial derivatives. In order to obtain a canonical Lagrangian, 
function of v{r,t), A(r, i) and their first derivatives, we use the method of Lagrange multipliers (l7L l29l |3Q| . Let us 
introduce the variable £ which is the Lagrange multiplier: this allows to set the Ricci scalar R as a constraint for the 
dynamics. In this way the point-like action can be written as 



S = / drdt r z e— i f(R) - £ 



R + - C 



1 2 2 ' 

v rr + -v r {v r — A r ) + -(y r — A r ) + -it 



-\ t (v t - A t ) - A M 



(4) 

By this approach, we can consider i?(r, t), z/(r, t) and A(r, t) as independent Lagrangian variables. The Lagrange 
multiplier can be easily found by varying the action with respect to R, and the result is 



Z = f'(R), 



(5) 



where the prime denotes the derivative with respect to R. Using this expression for £ and integrating by parts, action 
can be rewritten as 



S = 



dtdr r 2 e 



f(R) - Rf'(R) - 2 



f'(R) 



—A 



'f"{R)Rt\t 



(G) 



Clearly, action ([6]), when varied with respect to R, v, and A, gives the same equations of ([T]). Furthermore, we are 
disregarding a divergence and this fact could affect the conserved current. Finally, we are going to consider internal 
(i.e. vertical) transformations so that it is not important to know the transformation rules of v and A with respect to 
spacetime transformations. The equations of motion are 



-f"(R)R r + lf"(R)X r Rr - f"\R)R 2 r - f"(R)Rrr + ~f\R)K - \f'(R) 



+ e* 



\R!'{R) \f{R) + ^f'(R) 



-e-»f"(R)R t \ t = 0, 



(7) 



for A and 
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--f"(R)v t R t + f"(R)R 2 + f"(R)R t 



\f(R) ~ -f(R>r - -f"(R)Rr - \f"{R)R r v r 
r z r r I 



\Rf'{R)~; \f{R) + \f{R) 
2 2 H 



0. 



(8) 



for v. The equation for R(r,t) is given by Eq. ([3]), that is the Euler constraint on the dynamics. Finally, we obtain a 
system of three partial differential Eqs. (([7]), © and ©) in the unknown functions v(r,t), X(r,t) and the constraint 
R(r,t), once f(R) is given. As a first check it is easy to see that Eqs. ([7]) and allow us to recover the well-known 
results of GR by replacing f(R) = R. 

It is straightforward to show that when curvature is constant over space and time (i? = Ro), an interesting solution 
can be found. Indeed subtracting Eqs. (UJ and we immediately obtain 



A r + v r = 0, 

and substituting into Eq. ([7]), the solution for the metric component g rr is 



(9) 



e A (^) = - (10) 

1 3 1r 



where is an integration function depending on time, and 



A 



Ro 



f(Ro) ~ 
f'(Ro). 



(ii) 



Solving Eq. ([9]), we get 



= A 2 {t)e- x ^\ (12) 

where ^(t) is another integration function depending on time. Inserting Solutions (|10[) and (|12|) into Eq. ([3]), we 
have a constraint on the integration functions: it is easy to show that the only way to make Eq. ([3|) consistent is to 
choose A\{t) — A\ ^constant, for which we recover the well-known static Schwarzschild-De Sitter solution 



s tf < r '*> = A 2 (t) 



1 - 



T 



A 1 (t) 



2r 



ds 2 = 



Ax 
2r 



dP 



Ar 2 
3 



2r 



■dr l 



r 2 dQ?, 



where time as been scaled accordingly 



dP = A 2 (t)dP, (13) 

being Rq = — 4A. The metric is time independent, and we have thus proved the Birkhoff theorem in the case of 
/(i?)-gravity with constant curvature. It is worth noticing that such a theorem does not hold in general for any 
/(i?)-model as shown in (l3l. [l9|. 



III. NOETHER SYMMETRY APPROACH 



Let us consider now the Noether Symmetry Approach by which it is possible to obtain conserved quantities asking 
for the invariance of the Lagrangian under the Lie derivative along an appropriate vector field. In this way, constraints 
on dynamics are achieved and we it is possible to solve the equations of motion. In general, a point-like canonical 
Lagrangian C depends on the variables q J (x fJ ') and on their derivatives d v q 3 (x^). The corresponding Euler- Lagrange 
equations are 



4 



_dc_ _ dc 



Contracting eq. ()14ll with some unknown functions a? = a J (q l ), it gives 

DC dC 

Since we can write 



dd^qi dq 



<* J d.^TZT-inj =0- (15) 



from Eq. (fT5|) we immediately obtain 



a 'if^)^ d 4 + {d ^^ = L ^ (17) 

where Lx denotes the Lie derivative along the vector field 

which represent the generator of symmetry. We can immediately infer the Noether Theorem which states that if 
Lx£ = 0, the Lagrangian C is invariant along the vector field X. As a consequence, we can define the current [3lJ 

f = ai^f—, (19) 

which is conserved being 



d^f - 0. (20) 

As already developed in [34j . the presence of Noether symmetries allows to reduce dynamics and then to find out 
exact solutions. Specifically, symmetries select also the analytic form of f(R). A similar approach has been developed 
for cosmological solutions |27| . 



IV. NOETHER SYMMETRIES IN SPHERICALLY SYMMETRIC f{R) GRAVITY 

Let us take into account the case of spherical-symmetry in /(i?)-gravity. For the sake of generality, we can assume 
that the vector field X depends on the whole set of configuration variables that are functions of the radial and time 
coordinates. The configuration space is Q = {i/(r,t), X(r,t), R(r,t)}, so that the generator of symmetry 1 becomes: 

d „ d d d „ d d d „ d d 

X = a— +/3— + 7 — + a r — + 7r — +a t — +/3 t — + 7t — , 21 
ov o\ oR ov T o\ r oR r ov t a\ t oR t 

where a, /3 and 7 are functions of the variables v, X and R. As discussed above, a symmetry exists if the equation 
LyiC — has solutions for the functions a, (3 and 7 where at least one of them is different from zero. To verify such 



It is worth noticing that R, due to the constraint Eq. |3jl, is not a proper variable. However, it is needed in order to make canonical the 
Lagrangian. 
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a condition, we can set to zero the coefficients of terms R 2 , > v rX Tl v r R r , \ r R r , v t \t, vtRt and XtRt in 
Lx£ = [33. We obtain the following system of partial differential equations, linear in a, f3 and 7: 



rV^/"(i*)^=0, (22) 



r 2 e^/"(i*)||=0, (23) 



r 2 e^/"(i?)^ = 0, (24) 



r 2 e^/"(i?)||=0, (25) 



r 2 e^/"(i?)^ = 0, (26) 



a > da V+SS7 + |S = 0, (27) 



2 2" /"(i?) ' ai? 

-r 2 e^/"(i?)^ = 0, (28) 

r 2 e^/"(i?)^=0, (29) 

r 2 e^/"(i?)^=0, (30) 

2 e ^- a) + ^ + 7 7W + 0a ' (31) 

where the following constraint has to be satisfied 

ir 2 (f(R) - Rf'(R)) (a + /3)- f(R) (a + (3) - r 2 i?/"(i?) 7 - 2/"(i?) 7 + e~ A [/'(i?) (a - p) + 

-rX r f(R) (a - p) - 2rX r f(R)£- + 2 1 f"{R) (1 - r\ r )] = 0. (32) 

If we consider that f"(R) ^ 0, that is the trivial solution f(R) = R is excluded, Eqs. ([2"2 > ([2l) f and ([2"g ]) -([3"0" l) can be 
immediately solved giving the result 

a = a(y), /3 = ^(A) , 7 = -y(R) . (33) 

Eq. (|2T|) can be recast in the form 

a da 1 f"'(R) d 7 m . 
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where fj,o is a constant, because the r.h.s. of the first equality depends only on R while the l.h.s. one is a function of 
v and A, which are independent variables. We can assume /i = 0, so we obtain, from Eq. (|34|) the following set of 
equations: 



a da 1 
2 + d^ = 2 ft 



(35) 



(h, 

7 



f"{R) 



dR. 



(36) 



By the same arguments used above, we obtain, from Eq. (|35]), that j3 is constant, and using this result and putting 
Eq. fl2| into Eq. (J3TJ) , we get 

a = f3 = a , (37) 
where ao is a constant, being this result consistent with Eq. (|3"5")) . Moreover, integrating Eq. (|3"6"]) we obtain 
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Mi 



being /ii an integration constant. Now the constraint (|32[) becomes 



(38) 



a (/(#) - i?/'(i?)) r 2 - 2a f(R) - ^r 2 R - 2^1 + 2^^ (1 - r\ r ) = 0. 
Using Eq. (|20l) for the conservation of the current j^ 1 , we immediately obtain 



(39) 



(r 2 e^) U /"(i?)i?, +«o-/'(i?) + Mi^j = d t [(^V) (a f"(R)Rt + »i\t) 



For the particular choice ao = 0, the constraint ([39]) becomes 



(40) 



i? 



while Eq. (gOl) is 



(41) 



8 r 



(42) 



Using Eq. (|42p into the expression for the Ricci scalar (J3]), it is immediate to verify that the dynamical constraint 
(|4ip is satisfied. In this way we have that our system of equations is consistent having the four equations ([7]), (J5J, 
(|4Tj) and (|42|) . in the four unknown functions v(r,t), X(r,t), R(r,t) and f(R). Our task is now to determine these 
unknown functions selected by the presence of the Noether symmetry. 



V. EXACT SOLUTIONS 



We are able to solve the dynamics by the Noether symmetry approach. Considering Eqs. ([7]), © and (|4"Tj) . in the 
case of constant curvature R = Rq, we obtain the following solutions for the metric coefficients 



(43) 
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where Bi(t) and B 2 (t) are integration functions. Moreover we get as constraint f(Ro) — 0. Then Eq. (|4"2")l for the 
conserved current, using the solution (|4"3")) with Bi(t) — B^constant, gives 

R Q = 0. (44) 



Finally, using this result in Eq. (|43|) and rescaling the time as in Eq. (|T3|) . we obtain the following time-independent 
solution 

ds2 = ( x ~ dp jr^) dr2 ~ rW ' (45) 

which is the well-known Schwarzschild solution in GR. We can conclude that Noether Symmetry Approach for f(R)- 
gravity reproduces the same solutions of GR in the case of constant curvature. This is not a new result but it is useful 
to test the method [l7j . 

However, other interesting solutions can be found in the case of non-constant curvature in a spherically symmetric 
static metric. In fact, the current conservation Eq. (|42l) becomes now 



d / 



Ve^z^J = 0, (46) 
where v — v(r) and A = A(r). Assuming, as standard, that v = —A [33j, Eq. (|4"6"]) can be easily solved to give 

v{r) = In (Ct - -) + C 2 = -A(r), (47) 



r / 

being Cj and C*2 integration constants. Subtracting Eqs. ([7]) and ([5]), after some algebra, we obtain 

f'(R(r))=D 1 r + D 2 , (48) 

where, again, D\ and D 2 are integration constants. In order to find the form of the /(i?)-model, we need to express 
f'(R) as an explicit function of R. Then using the expression for the Ricci scalar (|41|) . which now is 

„„ . 2<^i) , (49, 
we obtain r as a function of R. Substituting into Eq. (|48p . we are able to perform the integration obtaining 

f(R) = 2Di (2C ie - C2 -2y R^ + D 2 R + D 3 , (50) 

where D3 is a constant. It can be demonstrated that our system of equations, considering also Eq. ©, is solved for 
given choices of the constants. For C\ = h, C 2 = 0, D 2 = — 3Di and D 3 = 0, we get the metric 

ds 2 = (l--)dt 2 - n 1 ^ dr 2 -r 2 dn 2 . (51) 

The selected /(i?)-function is 



2 rj (i-I) 



f(R)=D 2 R--D 



-R. 



(52) 



The Ricci scalar is uniquely determined as 



This solution is physically consistent for r > and it is asymptotically flat. 



(53) 
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VI. CONCLUSIONS 



In the context of /(-R)-gravity, a Lagrangian approach has been developed to study dynamics of spherically sym- 
metric metrics. We have obtained the Euler-Lagrange equations and solved them in the case of constant curvature 
R = R recovering the standard Schwarzschild-de Sitter solution of GR. In this case, the Birkhoff theorem holds. 
Then the Noether Symmetry Approach has been developed in order to reduce dynamics by finding out conserved 
quantities, which can be expressed as currents. After, we have solved the Noether system related to the condition 
Lx£ = 0, deriving a constraint equation over the dynamics and a conserved current expressed as functions of the 
dynamical variables. The constraint allows to select the form of the /(i?)-model. From this point of view the Noether 
symmetry ensures the closure of the dynamical system. Finally, we have presented some particular cases showing 
that Noether symmetries are compatible with GR and that solutions exist also for other /(i?)-models different from 
f(R) = R. In conclusion, the method reveals a useful approach both to select consistent /(-R)-models and to find out 
exact solutions. 
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